Physicochemical analysis of rotavirus segment 11 supports a 'modified panhandle' structure and not the predicted alternative tRNA-like structure (TRLS)

.Rotaviruses are a major cause of acute gastroenteritis, which is often fatal in infants. The viral genome consists of 11 double-stranded RNA segments, but little is known about their cis-acting sequences and structural elements. Covariation studies and phylogenetic analysis exploring the potential structure of RNA11 of rotaviruses suggested that, besides the previously predicted "modified panhandle" structure, the 5' and 3' termini of one of the isoforms of the bovine rotavirus UKtc strain may interact to form a tRNA-like structure (TRLS). Such TRLSs have been identified in RNAs of plant viruses, where they are important for enhancing replication and packaging. However, using tRNA mimicry assays (in vitro aminoacylation and 3'- adenylation), we found no biochemical evidence for tRNA-like functions of RNA11. Capping, synthetic 3' adenylation and manipulation of divalent cation concentrations did not change this finding. NMR studies on a 5'- and 3'-deletion construct of RNA11 containing the putative intra-strand complementary sequences supported a predominant panhandle structure and did not conform to a cloverleaf fold despite the strong evidence for a predicted structure in this conserved region of the viral RNA. Additional viral or cellular factors may be needed to stabilise it into a form with tRNA-like properties

Assessing the Utility of Thermodynamic Features for microRNA Target Prediction under Relaxed Seed and No Conservation Requirements

BACKGROUND: Many computational microRNA target prediction tools are focused on several key features, including complementarity to 5'seed of miRNAs and evolutionary conservation. While these features allow for successful target identification, not all miRNA target sites are conserved and adhere to canonical seed complementarity. Several studies have propagated the use of energy features of mRNA:miRNA duplexes as an alternative feature. However, different independent evaluations reported conflicting results on the reliability of energy-based predictions. Here, we reassess the usefulness of energy features for mammalian target prediction, aiming to relax or eliminate the need for perfect seed matches and conservation requirement. METHODOLOGY/PRINCIPAL FINDINGS: We detect significant differences of energy features at experimentally supported human miRNA target sites and at genome-wide sites of AGO protein interaction. This trend is confirmed on datasets that assay the effect of miRNAs on mRNA and protein expression changes, and a simple linear regression model leads to significant correlation of predicted versus observed expression change. Compared to 6-mer seed matches as baseline, application of our energy-based model leads to ∼3-5-fold enrichment on highly down-regulated targets, and allows for prediction of strictly imperfect targets with enrichment above baseline. CONCLUSIONS/SIGNIFICANCE: In conclusion, our results indicate significant promise for energy-based miRNA target prediction that includes a broader range of targets without having to use conservation or impose stringent seed match rules

RNA secondary structure prediction from multi-aligned sequences

It has been well accepted that the RNA secondary structures of most functional non-coding RNAs (ncRNAs) are closely related to their functions and are conserved during evolution. Hence, prediction of conserved secondary structures from evolutionarily related sequences is one important task in RNA bioinformatics; the methods are useful not only to further functional analyses of ncRNAs but also to improve the accuracy of secondary structure predictions and to find novel functional RNAs from the genome. In this review, I focus on common secondary structure prediction from a given aligned RNA sequence, in which one secondary structure whose length is equal to that of the input alignment is predicted. I systematically review and classify existing tools and algorithms for the problem, by utilizing the information employed in the tools and by adopting a unified viewpoint based on maximum expected gain (MEG) estimators. I believe that this classification will allow a deeper understanding of each tool and provide users with useful information for selecting tools for common secondary structure predictions.Comment: A preprint of an invited review manuscript that will be published in a chapter of the book `Methods in Molecular Biology'. Note that this version of the manuscript may differ from the published versio

MicroTar: predicting microRNA targets from RNA duplexes

BACKGROUND: The accurate prediction of a comprehensive set of messenger RNAs (targets) regulated by animal microRNAs (miRNAs) remains an open problem. In particular, the prediction of targets that do not possess evolutionarily conserved complementarity to their miRNA regulators is not adequately addressed by current tools. RESULTS: We have developed MicroTar, an animal miRNA target prediction tool based on miRNA-target complementarity and thermodynamic data. The algorithm uses predicted free energies of unbound mRNA and putative mRNA-miRNA heterodimers, implicitly addressing the accessibility of the mRNA 3' untranslated region. MicroTar does not rely on evolutionary conservation to discern functional targets, and is able to predict both conserved and non-conserved targets. MicroTar source code and predictions are accessible at , where both serial and parallel versions of the program can be downloaded under an open-source licence. CONCLUSION: MicroTar achieves better sensitivity than previously reported predictions when tested on three distinct datasets of experimentally-verified miRNA-target interactions in C. elegans, Drosophila, and mouse

Molecular Analysis of Virulent Determinants of Enterovirus 71

Enterovirus 71 (EV71) is the most important causative agent of hand, foot and mouth disease (HFMD) in children. In most cases, it is a self-limiting illness. However some EV71 infectious cases can develop severe clinical outcomes, such as encephalitis, meningitis, poliomyelitis like paralysis, and even death. To identify the determinants of virulence, the deduced amino acid sequence of polyprotein and nucleotide sequence of 5′-NTR and 3′-NTR in 25 SC-EV71 strains (strains from severe cases) and 31 MC-EV71 strains (strains from mild cases) were analyzed. Results showed four amino acids on two positions (GlyP710/GlnP710/ArgP710 and GluP729) on the DE and EF loop of VP1, one (LysP930) on the surface of protease 2A and four nucleotides on three positions (GP272, UP488 and AP700/UP700) in the 5'-NTR region are associated with EV71 virulent phenotype. Predicted secondary structure of RNA using the consensus sequence of 5'-NTR by RNAStructure showed the mutation of nucleotide at position 488 in strain BJ08-Z004-3 (position 491 in prototype strain BrCr) can result in the discrepancy of an additional pair of nucleotides and thus change the stability of the second structure of IRES. Fragment base content analysis showed that in the region 696 to 714 bp at the 5'-NTR, where the AP700/UP700 was located, the nucleotide constitution ratios differed significantly between SC-EV71 and MC-EV71 strains. In conclusion, comparative genomic analysis showed that virulence of EV71 strains are mainly determined by the amino acids on two positions of VP1, one position of protease 2A and the nucleotides on three positions in 5'-NTR

ViennaRNA Package 2.0

<p>Abstract</p> <p>Background</p> <p>Secondary structure forms an important intermediate level of description of nucleic acids that encapsulates the dominating part of the folding energy, is often well conserved in evolution, and is routinely used as a basis to explain experimental findings. Based on carefully measured thermodynamic parameters, exact dynamic programming algorithms can be used to compute ground states, base pairing probabilities, as well as thermodynamic properties.</p> <p>Results</p> <p>The <monospace>ViennaRNA</monospace> Package has been a widely used compilation of RNA secondary structure related computer programs for nearly two decades. Major changes in the structure of the standard energy model, the <it>Turner 2004 </it>parameters, the pervasive use of multi-core CPUs, and an increasing number of algorithmic variants prompted a major technical overhaul of both the underlying <monospace>RNAlib</monospace> and the interactive user programs. New features include an expanded repertoire of tools to assess RNA-RNA interactions and restricted ensembles of structures, additional output information such as <it>centroid </it>structures and <it>maximum expected accuracy </it>structures derived from base pairing probabilities, or <it>z</it>-<it>scores </it>for locally stable secondary structures, and support for input in <monospace>fasta</monospace> format. Updates were implemented without compromising the computational efficiency of the core algorithms and ensuring compatibility with earlier versions.</p> <p>Conclusions</p> <p>The <monospace>ViennaRNA Package 2.0</monospace>, supporting concurrent computations <monospace>via OpenMP</monospace>, can be downloaded from <url>http://www.tbi.univie.ac.at/RNA</url>.</p

Local Absence of Secondary Structure Permits Translation of mRNAs that Lack Ribosome-Binding Sites

The initiation of translation is a fundamental and highly regulated process in gene expression. Translation initiation in prokaryotic systems usually requires interaction between the ribosome and an mRNA sequence upstream of the initiation codon, the so-called ribosome-binding site (Shine-Dalgarno sequence). However, a large number of genes do not possess Shine-Dalgarno sequences, and it is unknown how start codon recognition occurs in these mRNAs. We have performed genome-wide searches in various groups of prokaryotes in order to identify sequence elements and/or RNA secondary structural motifs that could mediate translation initiation in mRNAs lacking Shine-Dalgarno sequences. We find that mRNAs without a Shine-Dalgarno sequence are generally less structured in their translation initiation region and show a minimum of mRNA folding at the start codon. Using reporter gene constructs in bacteria, we also provide experimental support for local RNA unfoldedness determining start codon recognition in Shine-Dalgarno–independent translation. Consistent with this, we show that AUG start codons reside in single-stranded regions, whereas internal AUG codons are usually in structured regions of the mRNA. Taken together, our bioinformatics analyses and experimental data suggest that local absence of RNA secondary structure is necessary and sufficient to initiate Shine-Dalgarno–independent translation. Thus, our results provide a plausible mechanism for how the correct translation initiation site is recognized in the absence of a ribosome-binding site

RNA targeting with CRISPR–Cas13

RNA has important and diverse roles in biology, but molecular tools to manipulate and measure it are limited. For example, RNA interference1-3 can efficiently knockdown RNAs, but it is prone to off-target effects4, and visualizing RNAs typically relies on the introduction of exogenous tags5. Here we demonstrate that the class 2 type VI6,7 RNA-guided RNA-targeting CRISPR-Cas effector Cas13a8(previously known as C2c2) can be engineered for mammalian cell RNA knockdown and binding. After initial screening of 15 orthologues, we identified Cas13a from Leptotrichia wadei (LwaCas13a) as the most effective in an interference assay in Escherichia coli. LwaCas13a can be heterologously expressed in mammalian and plant cells for targeted knockdown of either reporter or endogenous transcripts with comparable levels of knockdown as RNA interference and improved specificity. Catalytically inactive LwaCas13a maintains targeted RNA binding activity, which we leveraged for programmable tracking of transcripts in live cells. Our results establish CRISPR-Cas13a as a flexible platform for studying RNA in mammalian cells and therapeutic development.National Institute of Mental Health (U.S.) (Grant 5DP1-MH100706)National Institute of Mental Health (U.S.) (Grant 1R01-MH110049

Optimal Use of Conservation and Accessibility Filters in MicroRNA Target Prediction

It is generally accepted that filtering microRNA (miRNA) target predictions by conservation or by accessibility can reduce the false discovery rate. However, these two strategies are usually not exploited in a combined and flexible manner. Here, we introduce PACCMIT, a flexible method that filters miRNA binding sites by their conservation, accessibility, or both. The improvement in performance obtained with each of these three filters is demonstrated on the prediction of targets for both i) highly and ii) weakly conserved miRNAs, i.e., in two scenarios in which the miRNA-target interactions are subjected to different evolutionary pressures. We show that in the first scenario conservation is a better filter than accessibility (as both sensitivity and precision are higher among the top predictions) and that the combined filter improves performance of PACCMIT even further. In the second scenario, on the other hand, the accessibility filter performs better than both the conservation and combined filters, suggesting that the site conservation is not equally effective in rejecting false positive predictions for all miRNAs. Regarding the quality of the ranking criterion proposed by Robins and Press and used in PACCMIT, it is shown that top ranking interactions correspond to more downregulated proteins than do the lower ranking interactions. Comparison with several other target prediction algorithms shows that the ranking of predictions provided by PACCMIT is at least as good as the ranking generated by other conservation-based methods and considerably better than the energy-based ranking used in other accessibility-based methods